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Abstract 

The physics of polarized nucleon-nucleon collisions originating from an internal polar- 
ized target in the HERA proton beam is investigated. Based on 240 pb _1 integrated lumi- 
nosity at 40 GeV cm. energy, statistical sensitivities are given over a wide (xi?,pT)~range 
for a variety of inclusive and exclusive final states. By measuring single spin asymmetries 
unique information can be obtained on higher twist contributions and their p^-dependence. 
From double spin asymmetries in both photon and J/ip production it appears possible to 
measure the polarized gluon distribution in the range 0.1 < x g i uon < 0.4 with a good 
statistical accuracy. 
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1 Introduction 



There is a widespread general consensus, much grown in the last years and based both on 
surprising experimental results and intense theoretical activity, that the spin structure and the 
spin dynamical properties of nucleons and hadrons in general are far from being understood; 
a satisfactory knowledge of hadronic structure and dynamics and their correct description in 
terms of constituents cannot ignore the spin subtleties and more experimental information is 
badly required. 

An experiment ( i HERA-N y utilising an internal polarized nucleon target in the 820 GeV 
HERA proton beam would constitute a natural extension of the studies of the nucleon spin 
structure in progress at DES Y with the HERMES experiment . Conceivably, this would be 
the only place where to study high energy nucleon-nucleon spin physics besides the dedicated 
RHIC spin program at BNL || supposed to start early in the next decade. 
An internal polarized nucleon target offering unique features such as polarization above 80% 
and no or small dilution, can be safely operated in a proton ring at high densities up to 10 14 
atoms/cm 2 |3[]. As long as the polarized target is used in conjunction with an unpolarized 
proton beam, the physics scope of HERA-N would be focused to 'phase I', i.e. measurements of 
single spin asymmetries. Once later polarized protons should become available, the same set-up 
would be readily available to measure a variety of double spin asymmetries. These 'phase IF 
measurements would constitute an alternative - fixed target - approach to similar physics which 
will be accessible to the collider experiments STAR and PHENIX at the low end of the RHIC 
energy scale (y^ ~ 50 GeV) ||. 

We shall briefly discuss here the physics motivations to measure single and double spin asym- 
metries in several pN inclusive and exclusive processes; more details and further discussions, 
together with a complete list of references, which we cannot give here for lack of space, can be 
found in Ref. ||. We recall that the integrated luminosity calculation is based upon realistic 
figures. For the average beam and target polarisation Pb = 0.6 and Pt = 0.8 are assumed, 
respectively. A combined trigger and reconstruction efficiency of C ~ 50% is anticipated. Using 
Ib = 80 niA = 0.5 • 10 18 s" 1 for the average HERA proton beam current (50% of the design 
value) and a rather conservative polarized target density of tit = 3 • 10 13 atoms/cm 2 the pro- 
jected integrated luminosity becomes £ ■ T = 240 pb~ x when for the total running time T an 
equivalent of T = 1.6 • 10 7 s at 100 % efficiency is assumed. This corresponds to about 3 real 
years under present HERA conditions. However, proton currents much higher than the original 
HERA design value of 160 mA are envisaged in the HERA luminosity upgrade program. In 
addition, experience from UA6 running at CERN shows that after having gained some practical 
running experience it presumably becomes feasible to operate the polarized gas target at about 
3 times higher density without seriously affecting the proton beam lifetime. Hence in a few years 
500 pb -1 per year will presumably become a realistic figure. 

We note that, except in the case of single spin asymmetries, we took into account major accep- 
tance limitations and jet detection efficiencies. Hence it can be anticipated that the sensitivities 
shown in the rest of the paper are realistic for an about one year's running of a future polar- 
ized nucleon-nucleon scattering experiment at HERA. Any better estimate requires considerably 
intensified efforts to be invested along many different directions, like machine and target limita- 
tions, detector capabilities versus rate, acceptance, costs etc. 
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2 Single Spin Asymmetries 



Single spin asymmetries in large pt inclusive production, both in proton-nucleon and lepton- 
nucleon interactions, have recently received much attention (for references see ||). The naive 
expectation that they should be zero in perturbative QCD has been proven to be false, both 
experimentally and theoretically. It is now clear that higher twist effects are responsible for these 
asymmetries, which should be zero only in leading twist-2 perturbative QCD. 
Several models and theoretical analyses suggest possible higher twist effects: there might be 
twist-3 dynamical contributions, which we shall denote as hard scattering higher twists; there 
might also be intrinsic k± effects, both in the quark fragmentation process and in the quark 
distribution functions. The latter are not by themselves higher twist contributions - they are 
rather non-perturbative universal nucleon properties - but give rise to twist-3 contributions when 
convoluted with the hard scattering cross sections. The dynamical contributions result from a 
short distance part calculable in perturbative QCD with slightly modified Feynman rules, com- 
bined with a long distance part related to quark-gluon correlations. 

An intrinsic k± effect in the quark fragmentation is known as Collins or sheared jet effect; it sim- 
ply amounts to say that the number of hadrons h (say, pions) resulting from the fragmentation of 
a transversely polarized quark, with longitudinal momentum fraction z and transverse momen- 
tum k±, depends on the quark spin orientation. That is, one expects the quark fragmentation 
analysing power A q (k±) to be different from zero: 

where, by parity invariance, the quark spin should be orthogonal to the q — h plane. Notice also 
that time reversal invariance does not forbid such quantity to be 7^ because of the (necessary) 
soft interactions of the fragmenting quark with external strong fields, i.e. because of final state 
interactions. This idea has been applied to the computation of the single spin asymmetries ob- 
served in pp^ — > ixX 0. 

A similar idea applies to the distribution functions, provided soft gluon interactions between 
initial state partons are present and taken into account, which most certainly is the case for 
hadron-hadron interactions. That is, one can expect that the number of quarks with longitudinal 
momentum fraction x and transverse intrinsic motion k± depends on the transerve spin direction 
of the parent nucleon, so that the quark distribution analysing power N q (kj_) can be different 
from zero: 

N^^ i^fl^-^f^^ O (2) 

Jq/m {x, k ± ) + f q/N l {X, fcjj 

This effect also has been used to explain the single spin asymmetries observed in pp^ — > nX || . 
Note that both A q (k±) and N q (k±) are leading twist quantities which, when convoluted with 
the elementary cross- sect ions and integrated over kj_, give twist-3 contributions to the single 
spin asymmetries. 
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Each of the above mechanisms might be present and might be important in understanding twist- 
3 contributions; in particular the quark fragmentation or distribution analysing powers look like 
new non-perturbative universal quantities, crucial in clarifying the spin structure of nucleons. It 
is then of great importance to study possible ways of disentangling these different contributions 
in order to be able to assess the importance of each of them. We propose here to measure the 
single spin asymmetry 

d(J AB^CX _ d(J AB^CX 
da AB^CX + da ABl-,CX 



in several different processes AB 1 * — > CX which should allow to fulfil such a task. To obtain a 
complete picture we need to consider nucleon-nucleon interactions together with other processes, 
like lepton-nucleon scattering, which might add valuable information. For each of them we 
discuss the possible sources of higher twist contributions, distinguishing, according to the above 
discussion, between those originating from the hard scattering and those originating either from 
the quark fragmentation or distribution analysing power. 

• pN^ -> hX 

In this process all kinds of higher twist contributions may be present; this asymmetry alone 
could not help in evaluating the relative importance of the different terms. 

• pN^ -> 7X, piV T -> fi + i2~X, pNl -> jets + X 

Here there is no fragmentation process, and we remain with possible sources of non-zero single 
spin asymmetries in the hard scattering or the quark distribution analysing power. 

• INI -> hX 

In such a process the single spin asymmetry can originate either from hard scattering or from 
k± effects in the fragmentation function, but not in the distribution functions, as soft initial 
state interactions are suppressed by powers of a em . Moreover, this process allows, in principle, a 
direct measurement of the Collins effect, i.e. of the quark fragmentation analysing power, via a 
measurement of the leading-twist difference of cross sections for the production of two identical 
particles inside the same jet, with opposite k±. 

• INl -> 7X, 7 iV T -> 7 X, IN 1 ' -> fi+fi-X, INI -> jets + X 

Each of these processes yields a single spin asymmetry which cannot originate from distribution 
or fragmentation k± effects; it may only be due to higher-twist hard scattering effects, which 
would thus be isolated. 

It is clear from the above discussion that a careful and complete study of single spin asymmetries 
in several processes might be a unique way of understanding the origin and importance of higher 
twist contributions in inclusive hadronic interactions; not only, but it might also allow a deter- 
mination of fundamental non-perturbative properties of quarks inside polarized nucleons and of 
polarized quark fragmentations. Such properties should be of universal value and applicability 
and their knowledge might be as important as the knowledge of unpolarized distribution and 
fragmentation functions. 

— * 

In the following we discuss the capability of HERA-iV to investigate some of these processes. 
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Inclusive pion production p^p — > tt 0± X at 200 GeV exhibits surprisingly large single spin 
asymmetries, as it was measured a few years ago by the E704 Collaboration using a transversely 
polarized beam ||. For any kind of pions the asymmetry shows a considerable rise above 
xp > 0.3, i.e. in the fragmentation region of the polarized nucleon. It is positive for both tt + 
and 7T° mesons, while it has the opposite sign for ir~ mesons. The charged pion data taken in 
the 0.2 < pt < 2 GeV range were split into two samples at pr = 0.7 GeV/c; the observed rise 
is stronger for the high sample, as can be seen from fig. [I]. 
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Figure 1: Single spin asymmetry in inclusive pion production p^ + p — > 7r 0± + X measured 
by the E704 Collaboration J^/ and shown for two subregions ofpf. 



Contours characterizing different HERA-iV sensitivity levels (5An = 0.001, 0.01 and 0.05) for 
an asymmetry measurements in the reaction pp^ — > tt + X are shown in fig. Note that in the 
large pt region the contours calculated with big Ap^ x Ap? bins are appropriate, since usually 
a larger bin size is chosen where the statistics starts to decrease. We can conclude that the 
accessible pt values are significantly larger than those E704 had; the combined px dependence 
of all involved higher-twist effects can be measured with good accuracy (5A^ < 0.05) up to 
transverse momenta of about 10 GeV/c in the central region \xf\ < 0.2 and up to 6 GeV/c 
in the target fragmentation region. This corresponds to an almost one order of magnitude ex- 
tension in the px range in comparison to E704. The capability of HERA-iV to really prove a 
predicted pt dependence is shown in fig. |3|, where the curve was obtained assuming a non-zero 
quark distribution analysing power (cf. eq. §) according to Ref. ||. 
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Figure 2: Contours of the asymmetry sensitivity levels for n + production in the (pr, xf) plane. 
Lines of constant laboratory angles of the pion are shown. 



Inclusive direct photon production, pp^ — > jX, proceeds without fragmentation, i.e. the 
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Figure 3: Capability of HERA-N to discriminate predictions for different pr- 
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photon carries directly the information from the hard scattering process. Hence this process 
measures a combination of initial effects and hard scattering twist-3 processes. The first and 
only results up to now were obtained by E704 Collaboration flKj] showing an asymmetry com- 
patible with zero within large errors for 2.5 < pt < 3.1 GeV/c in the central region \xp\ ~ 0.15. 
The experimental sensitivity of HERA-iV was determined using PYTHIA 5.7 by simultaneous 
simulation of the two dominant hard subprocesses contributing to direct photon production, 
i.e. gluon-Compton scattering (qg — > 79) and quark-antiquark annihilation (qq — > jg), and 
of background photons that originate mainly from 7r° and r\ decays. It turns out that a good 
sensitivity (about 0.05) can be maintained up to px < 8 GeV/c. For increasing transverse mo- 
mentum the annihilation subprocess and the background photons are becoming less essential; we 
expect to be able to detect a clear dependence on pt, of the direct photon single spin asymmetry. 



Inclusive J/ip production was calculated in the framework of the colour singlet model 11 



Our calculations at HERA-iV energies || show an asymmetry less than 0.01 in the region 
\xp\ < 0.6, i.e. the effect is practically unobservable. 



3 Double Spin Asymmetries 

Perturbative QCD allows sizeable lowest order double spin asymmetries for various 2— >2 par- 
tonic subprocesses. Relying on the factorization theorem a rich spectrum of asymmetries at the 
hadronic level can be predicted which constitute the backbone of the RHIC spin physics pro- 
gram. When both incoming particles are longitudinally polarized, the insufficient knowledge of 
the polarized gluon distribution makes the predictions for double spin asymmetries All to some 
extent uncertain. Conversely, the measurement of All i n certain final states seems to be one of 
the most valuable tools to measure the polarized gluon distribution function in the nucleon. The 
presently most accurate way to do so is the study of those processes which can be calculated 
in the framework of perturbative QCD, i.e. for which the involved production cross sections 
and subprocess asymmetries can be predicted. Both direct photon (plus jet) and J/ip (plus jet) 
production are most suited because there are only small uncertainties due to fragmentation. 
In the following we discuss the corresponding capabilities of HERA-iV, operated in doubly po- 
larized mode ('Phase II'), to perform such measurements. 



Inclusive photon production with HERA-iV was the subject of a very recent study fT2|| . 
Basing on a NLO calculation rather firm predictions were obtained including an assessment of 
their theoretical uncertainties; the latter turned out to be of rather moderate size. In fig. |3] three 
different predictions for the asymmetry are shown in dependence on px and pseudorapidity 77, in 
conjunction with the attainable statistical uncertainty of HERA-iV. We note that the dashed 



line is rather close to the prediction of Ref. JT3], set A. As can be seen, there is sufficient sta- 



tistical accuracy up to transverse momenta of about 8 GeV/c to discriminate between different 
polarized gluon distribution functions (cf. fig f|a). At pr = 6 GeV/c there is sufficient accuracy 
to check the asymmetry prediction for photon pseudorapidities between -1.5 and 1.5 (cf. fig|]b). 
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Figure 4: Inclusive photon production: Double spin asymmetry vs. a) px and b) rj for the 
NLO 'valence set' of Ref. Hi]/ (full line), shown in conjunction with the HERA-N statistical 
sensitivity. The dotted line corresponds to set C of Ref. [LDj] and the dashed line is close to set 
A of Ref. 



In photon plus jet production the away-side jet is measured as well and the complete kine- 
matics of the 2— >2 subprocess can be reconstructed. In this case the asymmetry A LL can be 
directly related to the polarized gluon distribution if a certain subprocess can be selected. Using 
this approach photon plus jet production was discussed in Ref. as a tool to directly measure 
AG/G. The quark-antiquark annihilation subprocess is suppressed relatively to quark-gluon 
Compton scattering because of the lower density of antiquarks (of the polarized sea) compared 
to gluons (polarized gluons). The absolute statistical error of AG(x g )/G(x g ) was obtained as 



AG{x g ) = 5A LL 

[ G{x 9 ) J A DIS -a LL 1 ) 

Here Aois an d (111 are to be taken at the appropriate values of xf and x g , respectively. In 
calculating the r.h.s. of eq. [| we had to take into account the influence of the acceptance, as 
described in ||. 

In fig. |5] the calculated HERA-iV statistical sensitivity, on the present level of understanding, 
is shown vs. x g i uon in conjunction with predicted errors for STAR running at RHIC at 200 GeV 
cm. energy |I5|. The errors demonstrate clearly that in the region 0.1 < x g < 0.4 a significant 



result can be expected. This statement will very probably remain valid if once the systematic 
errors will have been estimated. As can be seen, the measurement of AG/G from photon plus 
jet production in doubly polarized nucleon-nucleon collisions at HERA can be presumably per- 
formed with an accuracy being about competitive to that predicted for RHIC. 
At this point we note that the HERA-iV fixed target kinematics causes additional problems for 
the jet reconstruction when compared to a collider experiment. As obtained from rather prelimi- 
nary investigations 0, the number of photon events accompanied by a successfully reconstructed 
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jet decreases considerably when approaching lower values of and, correspondingly, of x g i uon . 
To have a more realistic statistical significance than that shown in Ref. H we now include 
these preliminary jet reconstruction efficiencies. As a result, the statistical error bars became 
somewhat larger for smaller transverse momenta. 
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Figure 5: Typical predictions for the polarized gluon distribution confronted with the projected 
statistical errors expected for HERA-N and RHIC experiments. 



J /ip Production. Compared to direct photon production the production of quarkonium states 
below the open charm threshold is a similarly clean tool to measure the polarized gluon distribu- 
tion. Hence many statements made in the previous section apply here, as well, and the principle 
of analysis is very similar. Because of the relatively large quark mass the cc production cross 
section and the expected asymmetry are supposed to be calculable perturbatively. 
Quarkonium production has traditionally been calculated in the color-singlet model (CSM) [ |16H 
where the quark-antiquark pair is produced in a color-singlet state with the quantum numbers 
of the corresponding hadron. This heavy mass pair then creates the hadronic state with a prob- 
ability determined by the appropriate quarkonium wave function at the origin. It is assumed 
that for heavy quarks soft gluon emission is negligible, as also other non-perturbative effects like 
higher twist contributions. While this model gives a reasonable description of J/ip production 
cross section shapes over p T and xp, it completely fails in the explanation of the integrated cross 
section; a K factor of 7 -r- 10 is needed to explain the data. The anomalously large cross section 



n] for J / ip production at large transverse momenta found at the Tevatron reveals another bad 
feature of the CSM; it is not able to explain the large ip' and direct J/ip production rates at 
CDF. All these observations led to the understanding that fragmentation and hadronization of 



color-octet |L8| qq pairs are essential in the heavy quarkonium production process. 

Within the framework of the color-octet mechanism the quarkonium production process can be 

separated, according to the factorization hypothesis, into a short and a long distance part. The 
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former describes the quark- ant iquark pair production at small distances and can be computed 
perturbatively. The latter is responsible for the creation of a particular hadronic state from the 
quark-antiquark pair; its matrix elements can not be calculated perturbatively. 
The shapes of the distribution of short distance matrix elements calculated within the color- 
octet model indicate that the new mechanism seems to be able to explain the Tevatron data for 
direct J/ip and ip' production at large pt |19||. We note that unlike color-singlet matrix elements 



being connected to the subsequent hadronic non-relativistic wave functions at the origin, color- 
octet long distance matrix elements are unknown and have to be extracted from the experiment. 
As it was shown recently [f20], the color-octet mechanism gives the dominant contribution in 
J/V'-production at HERA-iV energies. This concerns not only the total cross section | 20|| , but 



also J/ij) production at non-zero pt, i.e. at those values of transverse momenta which are not 
due to internal motion of partons inside the colliding hadrons, pt ~ 1.5 GeV/c [PH| . 

For inclusive J/ip production we calculated the double spin-asymmetry within the frame- 



work of the color-octet model using the value of the color-octet matrix elements from It is 
interesting to note that the magnitude of the asymmetry does almost not depend on the choice 
of the color-octet matrix elements, if one assumes that the transition matrix elements of octet 
3 Pj-states into J/ip are negligible compared to those of the 1 S'o state. This is the most likely 
scenario since only in this case it is possible to establish a consistence between different combi- 

or 



nations of the above two matrix elements, as extracted from CDF data ]T9[ and photo- |22 
hadroproduction [^(J data. We underline that the measurement of the double spin asymmetry in 
J/ip production would allow to extract information about the color-octet matrix elements sepa- 
rately, whereas from unpolarized experiments it is only possible to extract combinations of them. 
In fig. |0a we present the expected asymmetry versus px for two different sets of polarized gluon 
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Figure 6: Inclusive J/ip production: Double spin asymmetry vs. a) px and b) rj for the LO set A 
(full line) and the set C (dotted line) of Ref. fcT$, shown in conjunction with the HERA-N 
statistical sensitivity. 



distributions taken from JTJJ; the solid curve corresponds to set A and the dashed one to set C. 
For set A the asymmetry appears sufficiently large to be observed and its measurement would 
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allow to extract information about the polarized gluon distribution function in the nucleon. As 
can be seen from fig. |6|b, a very good discrimination between set A and set C is possible over the 
whole HERA-iV pseudorapidity interval. For the mass of the charm quark we chose m c = 1.5 
GeV/c 2 , as it was used for extraction of the color-octet matrix elements from experimental data 
|L9| , |20| . We found that unlike the J/ip production cross section the asymmetry does practically 
not depend on the charm quark mass; if we vary the charm quark mass from 1.35 to 1.7 GeV/c 2 
the magnitude of the asymmetry changes by about 3 -f- 5% over the whole considered pt region. 



In J /ip plus jet production the study of the double spin asymmetry would allow to access di- 
rectly the polarized gluon distribution function, similar to the case of photon plus jet production. 
For the absolute statistical error of AG(x g )/G(x g ) an expression similar to eq. |] is obtained: 

AgCgJ 5A LL 

[ G(x g ) J [AG/ G] -a LL 1 1 

with ^ = AaQg _> j/i/j g ) 

v(j99 -> J/^ 9) + [q(x)/G(x)] ■ a(gq -> J/tp q)' 

Here the quark-gluon subprocess was neglected, since its contribution to cl^l amounts only to 
about 10% compared to the gluon-gluon fusion subprocess. 

Unfortunately, taking into account the acceptance limitations as it was done for photon plus jet 
production, the J/ if) production cross section for HERA- N decreases significantly. Following the 
same principle of analysis as in Ref. ||, it turns out that the measurement of AG(x)/G(x) in 
J/ijj plus jet production is feasible only for x g = 0.1 -4- 0.2, i.e. for J/ip transverse momenta 
of about 2.5 GeV/c. This prediction is shown as an additional entry in fig. |^. Although being 
a single point only, we underline that this is a very important measurement, because the lowest 
lying point from photon plus jet production is obtained for rather small values of px where 
perturbative QCD is not expected to give reliable predictions. We note that the nature of the 
gluon-gluon subprocess has rather similar consequences for jet plus jet production at RHIC; the 
prediction |TSJ consists of only one single point at a similar value of x g , as well (cf. fig. []). 



4 Elastic Scattering 

Large unexpected spin effects in singly polarized proton-proton elastic scattering p + p^ — > p + p 
have been discovered many years ago and remain unexplained up to now. The single spin 
asymmetry A jy was found significantly different from zero in the region 1 ~ pf, ~ 7 (GeV/c) 2 
as it is shown in fig. [?| in conjunction with the HERA-iV statistical errors. At HERA-iV energies 
the detection of the recoil proton for p\ values in the range 5 -4- 12 (GeV/c) 2 requires a very 
large angular acceptance (up to 40 degrees) The forward protons for the same interval in 
p\ have laboratory angles of the order of a few milliradians and require a dedicated detector 
very close to the beam pipe. Note that although the accessible p\ range would be similar to the 
range explored at low energies the cm. scattering angle detected amounts to a few degrees at 
HERA-iV energies, only. 

The transverse single-spin asymmetry in elastic pp scattering at HERA-iV and RHIC energies 



has been calculated in a dynamical model that leads to spin-dependent pomeron couplings [23[ . 
The predicted asymmetry is about 0.1 for p\ = 4-4-5 (GeV/c) 2 with an expected statistical 
error of 0.01 4- 0.02 for HERA-iV, i.e. a significant measurement of the asymmetry A^ can be 
performed to test the spin dependence of elastic pp scattering at high energies. 
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Figure 7: Compilation of experimental data on the asymmetry in elastic proton-proton scattering 
as a function of p\. In addition the projected statistical errors attainable with HERA-N are 
shown. 



5 Conclusions 

The physics potential of polarized nucleon-nucleon collisions originating from an internal target 
in the 820 GeV HERA proton beam has been investigated. Single spin asymmetries, accessible 
already with the existing unpolarized beam, are found to be an almost unique and powerful tool 
to study the nature and physical origin of twist-3 effects; even more so when taken in conjunction 
with results of other experiments at HERA. When measuring the polarized gluon distribution 
through double spin asymmetries in photon (plus jet) and J /if) (plus jet) production - requiring 
a polarized HERA proton beam - the projected statistical accuracies are found to be comparable 
to those predicted for the spin physics program at RHIC. In addition, significant results can be 
obtained on the long-standing unexplained spin asymmetries in elastic scattering. 
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